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Abstract 


Proteins are biological polymers composed of linear chains of 20 different amino acids. The sequence of amino 
acids for every protein is unique, and guides its folding into intricate 3-dimensional shapes, known as protein folds. 
The TIM barrel is one such fold, and is characterized by an interior 8-stranded B-barrel, surrounded and enclosed 
by 8 a-helices. TIM barrels are named after triose phosphate isomerase (TIM), an enzyme first structurally charac- 
terized in 1975, which lends its name to the fold. TIM barrels are prevalent in all forms of life, and across diverse 
metabolic pathways, with over 10% of all enzymes adopting this fold. The majority of TIM barrels are thought to 
have evolved from a common ancestor through gene duplication and domain fusion processes. TIM barrels have 
been created by protein engineers using preexisting half-barrel templates and de novo, without an existing tem- 
plate. This review will discuss the topological, structural, evolutionary, and design characteristics of TIM barrels in 


detail. 


Introduction 


The TIM barrel is a structurally conserved protein fold, 
named after triose phosphate isomerase, a constituent 
enzyme of the glycolysis pathway whose structure was 
first solved in 1975." TIM barrels contain 200-250 
amino acid residues,! which form 8 B-strands and 8 a- 
helices. The B-strands are arranged into a parallel B- 
barrel, and are surrounded by the 8 a-helices. The inner 
B-barrel is in many cases stabilized by intricate salt- 
bridge networks.! Loops at the C-terminal ends of the 
B-barrel are responsible for catalytic activity!"! while 
N-terminal loops are important for the stability of the 
TIM-barrels. Structural inserts ranging from extended 
loops to independent domains may be inserted in place 
of these loops or at the N/C-terminals. TIM barrels are 
ubiquitous, with approximately 10% of all enzymes 
adopting this fold.!© Further, 5 of 7 enzyme commission 
(EC) enzyme classes include TIM barrel proteins.!”7!*! 
The TIM barrel fold is evolutionarily ancient, with many 
of its members possessing little similarity today,'*! in- 
stead falling within the twilight zone of sequence simi- 
larity.2°0-4 TIM barrels appear to have evolved through 
gene duplication and domain fusion events of half-bar- 
rel proteins,"2! with a majority of TIM barrels originating 
from acommon ancestor. This lead many TIM barrels to 
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possess internal symmetries."*! Further gene duplica- 
tion events of this ancestral TIM barrel lead to diverging 
enzymes possessing the functional diversity observed 
today. TIM barrels have also been a longstanding target 
for protein designers. Successful TIM barrel designs in- 
clude both domain fusions of existing proteins and de 
novo designs. Domain fusions experiments have re- 
sulted in many successful designs, “*#J45IH647108)191120) 
whereas de novo designs only yielded successes after 28 
years of incremental development. 


Structure 
Topology 


The TIM barrel gets its name from the enzyme triose 
phosphate isomerase (TIM), which was the first protein 
possessing the fold to be crystallized. TIM barrels (Fig- 
ure 1A) contain 200-250 amino acid residues,”! folded 
into 8 a-helices and 8 B-strands. The B-strands are ar- 
ranged into a parallel B-barrel. The defining property of 
TIM B-barrels is that they always possess a shear num- 
ber of 8.! The shear number is determined by picking a 
residue x on B-strand-1, and moving along the B-barrel, 
in a perpendicular direction to the direction of the 
strands, until residue y on the original B-strand-1 is 
reached. The number of residues between the start and 
end positions (|y-x|) is the shear number.”! Since the 
number of strands is equal to the Shear number, side- 
chains point alternatively towards the pore and the 
core, giving a 4-fold symmetry. The a-helices surround 
and completely enclose the inner B-barrel. 
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Figure 2 | Triose phosphate isomerase (TIM) isolated from chicken muscles (PDB: 1TIM), the archetypal TIM barrel 
enzyme. (A) Cartoon representation of the TIM barrel structure. a-helices are colored teal, B-strands are colored orange, 
and loops are colored green. Note that the C-terminal ends of B-strands are depicted with arrowheads. (B) Core and pore 
regions are highlighted. Amino acid residues belonging to the pore are colored blue. Amino acid residues belonging to 
the core are colored orange. Note that the TIM barrel is depicted in a top-down view, where the C-terminal ends of the 
B-barrel are pointed towards the reader. 
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Figure 2 | TIM barrel topology. a-helices are colored teal, loops are colored green, and B-strands are colored in two 
shades of orange. Lighter shades indicate residues pointing inward, towards the barrel pore. Darker shades indicate 
residues pointing outward, towards the barrel core. Cyan lines depict an example backbone B-barrel hydrogen bonding 
network. Note that side-chain hydrogen bonding networks are not depicted here. Interior B-barrel residues (pore residues) 
display a 4-fold geometric symmetry, despite emerging from an 8-strand B-barrel. This symmetry is illustrated as two 
example "layers" in red and blue. Each layer contains 4 residues that point towards the pore, and lie on the same plane 
perpendicular to the barrel axis. The shear number for TIM barrels is always 8, and is illustrated in magenta. Some TIM 
barrels naturally adopt, or are designed to adopt, two or four-fold symmetry. Example asymmetric units are also high- 
lighted. This figure has been adapted with permission from previously published work.(22] 


2 of 10 | WikiJournal of Science 


36 


Short loops typically connect the a and B secondary 
structures, forming a (Ba)s repeat topology. In some 
cases, structures ranging from extended loops to inde- 
pendent domains may be inserted in place of these 
loops, or may be attached to the N/C-terminals. All TIM 
barrel enzymes possess catalytic sites at the C-terminal 
end of the B-barrel,'*! and structural inserts present 
close to this end may aid in catalytic activity. The overall 
topology of the TIM barrel is illustrated in Figure 2. 


Core and pore regions 


TIM barrels contain two distinct buried regions, where 
amino acid residues are completely enveloped by their 
neighbors and lack access to solvent. We have previ- 
ously referred to these regions as the ‘core’ and 
‘pore',”2] and are depicted in Figure 1B. It should be 
noted that the term ‘pore’ is a misnomer, as no solvent 
channels exist within this region. The core region con- 
sists of all residues constituting the a-B interface, and 
lies exterior to the central B-barrel. The pore region con- 
sists of all interior B-barrel residues, which are sur- 
rounded and enclosed by the B-barrel backbone. 


Due to the pleated nature of B-strands, alternate resi- 
dues along a strand are almost evenly split between the 
pore (53%) and core (47%).!7! For B-barrels, 95% of their 
core residues are buried. Only 11% of their core residues 
are polar, possessing an affinity for water, and pos- 
sessing the ability to form hydrogen bonds or salt 
bridges. Similarly, 84% of B-strand pore residues are 
buried. However, 42% of their pore residues are polar. 
These residues form intricate salt bridge networks to 
compensate for their lack of solvent accessibility. 


TIM barrel stabilizing elements 


Salt bridges within TIM barrel pores are thought to con- 
tribute to the overall stability of the fold. An example of 
a large salt bridge network can be found in 2-deoxyri- 
bose-5-phosphate aldolase (Figure 3). This network was 
found to be conserved across the Class | aldolase fam- 
ily.2) 


The exact reason for the overrepresentation of polar 
residues and salt bridges within the pore remains un- 
clear. One study proposes that they improve foldability 
rather than thermodynamic stability of TIM barrels.!27) 
During the folding process, inner pore residues on B- 
strands would be exposed to water. Partially-folded 
BaBa modules, called foldons, would be energetically 
stabilized by polar pore residues during this stage of 
folding. 
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Figure 3 | Example salt bridge network in 2-deoxyribose- 
5-phosphate aldolase (PDB: 1P1X). Interactions are shown 
as cyan dashed lines. Polar residues are colored green. 
Polar amino acids aspartate (D), glutamate (E), lysine (K), 
and arginine (R), are shown here. 


In another study involving the S. solfataricus indole-3- 
glycerol phosphate synthase TIM barrel protein, a con- 
served BaBaB module was found to be an essential fold- 
ing template, which guided the folding of other second- 
ary structures.'5! B-barrel closure only occurred at the 
end of the folding process. In this case however, the au- 
thors credited branched aliphatic amino acids (valine, 
leucine, and isoleucine) for foldon stability. 


Another stabilizing element in TIM barrels is the B-hair- 
pin clamp. Side chain H-bond donors at the N-termini 
of even-numbered B-strands often form H-bonds with 
main chain amide hydrogens in preceding odd-num- 
bered B-strands.!°271 These clamps (or hydrophobic 
side chain bridge analogs) are conserved in 3 indole-3- 
glycerolphosphate synthase TIM barrel orthologs from 
the bacterial and archaeal kingdoms,'*! implying they 
arose in their last common ancestor and have been pre- 
served for over a billion years. 


Structural inserts 


The N/C-terminal and loop regions on TIM barrel pro- 
teins are capable of hosting structural inserts ranging 
from simple secondary structural motifs to complete 
domains. These domains aid in substrate recognition 
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Figure 4 | Examples of structural inserts at TIM barrel loop and N/C-terminal regions. (A) The Bacillus subtilis 
Orotidine 5'-monophosphate decarboxylase (PDB: 1DBT). Orotidine 5'-monophosphate is colored green. a-helical inserts 
are colored teal. The catalytic arginine residue (R215) is displayed as sticks. (B) Mycobacterium tuberculosis bifunctional 
histidine/tryptophan biosynthesis isomerase (PriA) (PDB: 2Y85). CdRP, the product of the TrpF reaction, is colored green. 
B-strand/loop interchangeable structures are colored orange. (C) Lactococcus lactis dihydroorotate dehydrogenase A 
(DHODA) (PDB: 2DOR). B-strands forming a sheet are colored orange. Extended loops are colored green. The cavity 
formed by these structures is displayed as a blue mesh. The product orotate is colored magenta. the cofactor FMN is 
colored pink. (D) Methylophilus methylotrophus trimethylamine dehydrogenase (PDB: 2TMD).!791 The Rossmann fold 
domain is colored according to secondary structural elements. Cofactor FMN is colored magenta. The [4Fe-4S]* is 


colored red. Note that substrate/product were not crystallized. 


and catalytic activity. Four diverse examples of TIM bar- 
rels containing additional motifs and domains are dis- 
cussed below. 


Bacillus subtilis Orotidine 5'-phosphate decarbox- 
ylase®?0! (PDB: 1DBT, Figure 4A) is a TIM barrel protein 
displaying 4 a-helices in place of the Ba loops typically 
present at the C-terminal of the B-barrel (residues 35- 
42, 89-91, 126-133, and 215-219). One of these helices 
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(R215—K219) contains a conserved arginine residue 
(R215) required for interacting with a phosphate moiety 
on orotidine 5'-monophosphate. The other helices were 
not found to host residues critical for catalytic activity, 
and may serve in structural roles. 


Mycobacterium tuberculosis bifunctional _histi- 
dine/tryptophan biosynthesis isomerase (PriA) (PDB: 
2Y85, Figure 4B) possesses the ability to catalyse two 


fC) 


reactions: (i) HisA reaction: the conversion of N-[(5- 
phosphoribosyl) formimino]-5-aminoimidazole-4-car- 
boxamide ribonucleotide (ProFAR) to N-[(5-phos- 
phoribulosyl)formimino]-5-aminoimidazole-4-carbox- 
amide ribonucleotide (PRFAR), and (ii) TrpF reaction: 
N-(5'-phosphoribosyl)-anthranilate (PRA) to 1-(O-car- 
boxyphenylamino)- 1'-deoxyribulose-5'-phosphate 
(CdRP).24 PriA is a TIM barrel enzyme that accommo- 
dates both substrates using active site loops (loops 1, 5, 
and 6, extended Ba loops at the C-terminal end of the 
B-barrel) that change conformation depending on the 
reactant present.2] Loop 1 wraps over the active site 
only in the presence of ProFAR. Loop5 wraps over the 
active site, adopting a B-sheet conformation in the 
presence of CdRP, or a knot-like conformation in the 
presence of ProFAR. Loop 6 wraps over the active site 
for all reactants. 


Lactococcus lactis Dihydroorotate dehydrogenase A 
(DHODA) (PDB: 2DOR, Figure 4C) is an example of a 
TIM barrel possessing B-sheets and extended loops 
over the C-terminal end of the B-barrel. DHODA cata- 
lyzes the oxidation of dihydroorotate to orotate, 771%"! 
which is part of the de novo uridine 5'-monophosphate 
(UMP) synthesis pathway. This oxidation is mediated by 
flavin mononucleotide (FMN). Here, B-sheets and ex- 
tended loops enclose the active site forming a cavity, 
while also hosting several catalytic residues. 


The Methylophilus methylotrophus trimethylamine de- 
hydrogenase”®! (PDB: 2TMD, Figure 4D) TIM barrel is 
an example of a complete domain insertion. Here, a 
Rossmann fold domain is inserted at the C-terminal end 
of the TIM-barrel. Trimethylamine dehydrogenase cat- 
alyzes the conversion of trimethylamine to formalde- 
hyde.®*) This reaction requires both a reduced 6-S-cys- 
teinyl Flavin mononucleotide (FMN) cofactor and a re- 
duced iron-sulphur ([4Fe-4S]*) center. FMN is cova- 
lently bound within the C-terminal region of the B-bar- 
rel. The [4Fe-4S]* center is too large to be accommo- 
dated within the TIM barrel, and is instead placed in 
close proximity, 7 A away, at the interface between the 
TIM barrel and Rossmann fold domains. 


Folding mechanisms 


The conservation of the TIM barrel fold is mirrored by 
the conservation of its equilibrium and kinetic folding 
mechanisms in bacterial paralogs with phylogenetically 
distinct lineages. Chemical denaturation of several nat- 
ural?I371 and 2 designed TIM barrel variants®” invaria- 
bly involves a highly populated equilibrium intermedi- 
ate. The kinetic intermediates that appear after dilution 
from highly denaturing solutions involve an early mis- 
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Figure 5| The reaction coordinate diagram for SsIGPS at pH 
7.8 and 25°C. The refolding reaction begins in the unfolded, 
U state, initially misfolds to the Isp intermediate state, par- 
tially unfolds to reach the la intermediate state whose conver- 
sion to the subsequent lg intermediate state is rate-limiting. 
The final step is the conversion of Ip to the native state, N. The 
Ia and lg kinetic intermediates correspond to the intermediate 
observed in equilibrium unfolding studies. The ordinate rep- 
resents the free energy of each state in the folding reaction 
mechanism in kcal mol -+. The abscissa represents the de- 
pendence of the difference in free energy between 2 states on 
the denaturant concentration and is proportional to the 
change in buried surface, referenced to the U state. The ki- 
netic folding mechanism, illustrating the flow of the unfolded 
protein to the native conformation is shown beneath the re- 
action coordinate diagram. 


folded species that must at least partially unfold to ac- 
cess the productive folding pathway.?°!97! The rate-lim- 
iting step in folding is the closure of the 8-stranded B- 
barrel, with the preceding, open barrel form corre- 
sponding to the equilibrium intermediate(Figure 5). 
Native-centric molecular dynamics simulations recapit- 
ulate the experimental results and point the way to test- 
able computational models for complex folding mecha- 
nisms. 29! 


Conserved fitness landscapes 


TIM barrel proteins possess an unusually high sequence 
plasticity, forming large families of orthologous and pa- 
ralogous enzymes in widely divergent organisms. This 
plasticity suggests a sequence landscape that allows for 
protein adaptation to a variety of environmental condi- 
tions, largely independent of phylogenetic history, 
while maintaining function. Chan et al.’°! used a deep 
mutational scanning approach and a competition as- 
say'“°lto determine the fitness of all possible amino acid 
mutants across positions in 3 hyperthermophilic indole- 
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Figure 6 | Experimentally derived fitness landscapes mapped from point mutations represent single steps from WT se- 
quence. Despite significant divergence of WT in sequence space, the fitness landscapes of IGPS orthologues remain correlated 
(dashed lines). Rather than traditional two-dimensional heatmaps, fitness values are displayed on a three-dimensional pinwheel, 
highlighting the wide range of possible fitness effects of a single sequence step. The profiles of the pinwheels are similar, indi- 
cating the correlation of fitness landscapes, even if WT sequences (centers of the wheels) are only 40% identical and widely 
separated. Principal component analysis demonstrates a correlation between experimental fitness landscapes and amino-acid 


preferences in evolved sequences. 


3-glycerolphosphate synthase (IGPS) TIM barrel en- 
zymes in supporting the growth of a yeast host lacking 
IGPS. Although the 2 bacterial and 1 archaeal IGPS en- 
zymes were only 30-40% identical in sequence, their fit- 
ness landscapes were strongly correlated: the same 
amino acids at the same positions in the three different 
proteins had very similar fitness. The correlation can be 
thought of as the conservation of the fitness landscape 
for a TIM barrel enzyme across evolutionary time (Fig- 
ure 6). 


Evolution and origins 


The predominant theory for TIM barrel evolution in- 
volves gene duplication and fusion, starting with a half- 
barrel that eventually formed a full TIM barrel. Multiple 
studies support the theory of divergent evolution from 
a single ancestor, and are discussed below. 
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Evolution from a common ancestor 


In the early 1990s, Farber et al.) noted that all TIM 
barrel structures solved at the time were enzymes, indi- 
cating divergence from a common ancestor. Further, all 
TIM barrels possessed active sites at the C-terminal end 
of B-barrels. Brandén'*! suggested that a common 
phosphate binding site, formed by a small a-helix and 
TIM barrel loops-7/8, strongly indicated divergent evo- 
lution. Copley et al."“1! further studied these phosphate 
groups, concluding that 12 of 23 SCOP (structural clas- 
sification of proteins)! TIM barrel families diverged 
from a common ancestor. Similarly, Nagano et al.'°! 
concluded that there were hints for common ancestry 
for 17 of the 21 CATH (classification of protein domain 
structures)?! TIM barrel families. Based on these re- 
ports, it is plausible that the majority of TIM barrel pro- 
teins evolved from a common ancestor. 
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Origin through gene duplication and do- 
main fusion 


Many TIM barrel proteins possess 2-fold, 4-fold or 8- 
fold internal symmetry, suggesting that TIM barrels 
evolved from ancestral (Ba)s, (Ba)2z, or Ba motifs 
through gene duplication and domain fusion. A good 
example of 2-fold internal symmetry is observed in the 
enzymes ProFAR isomerase (HisA) and imidazole glyc- 
erol phosphate synthase (HisF) of the Thermotoga mar- 
itima histidine biosynthesis pathway."7! They catalyze 2 
successive reactions in the pathway, possess 25% se- 
quence homology, and possess root-mean-square devi- 
ations (RMSDs) between 1.5-2 A, suggesting diver- 
gence from a common ancestor. More interestingly, the 
loops on the C terminal ends of both HisA and HisF 
showed a twofold repeated pattern, suggesting that 
their common ancestor also possessed 2-fold internal 
symmetry. Using these observations, Lang et al." con- 
structed a model for the evolution of the TIM barrels 
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(Figure 5). An ancestral half-barrel would have under- 
gone a gene duplication and fusion event, resulting ina 
single protein containing two half-barrel domains. 
Structural adaptations would have occurred, resulting 
inthe merging of these domains to form a closed B-bar- 
rel, and forming an ancestral TIM barrel. Functional ad- 
aptations would have also occurred, resulting in the 
evolution of new catalytic activity at the C terminal end 
of the B-barrel. At this point, the common ancestor of 
HisA and HisF would have undergone a second gene du- 
plication event. Divergent evolution of the duplicated 
genes of the ancestral TIM barrel would have resulted in 
the formation of HisA and HisF. 


Interestingly, this evolutionary model has been experi- 
mentally validated using rational protein design and di- 
rected evolution. Hocker et al. first fused two C-termi- 
nal halves of HisF, yielding HisF-CC. This construct was 
then stabilized by the insertion of an internal salt- 


Ancestral TIM barrel 


Structural & 
functional 
adaptation 


Figure 7 | Model for the evolution of TIM barrels through gene duplication and domain fusion, as proposed by Lang et al." 
This model described the evolution of enzymes HisA and HisF of the histidine biosynthesis pathway. Two gene duplication steps 
are thought to have occurred. The first gene duplication resulted in two half-barrels that later fused and evolved into an ancestral 
TIM barrel. The second gene duplication event lead to diversification, and the evolution of different TIM barrel enzymes catalyz- 


ing different reactions. 
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bridge, yielding HisF-C*C."° Seitz et al.@*! and Hécker 
et al.">! then stepwise further stabilized and solubilized 
HisF-C*C by optimizing the half-barrel interface, gener- 
ating HisF-C**C and HisF-C***C, respectively. The 
crystal structure of HisF-C***C revealed a 2-fold sym- 
metric TIM barrel, validating the possibility of natural 
domain fusion. Moreover, Hécker created the first chi- 
meric HisAF and HisFA TIM barrels using HisA and HisF 
half-barrels."°! These experiments led to the proposal 
of a novel means of diversification and evolution of TIM- 
barrel enzymes through the exchange of (Ba)4 half-bar- 
rel domains amongst preexisting TIM barrels. In accord- 
ance with this idea, Claren et al. established high cata- 
lytic activity on the HisAF construct.” Similarly, chi- 
meric Bas-flavodoxin-like fold (CheY)/HisF TIM bar- 
rels,8129] and a perfectly 2-fold symmetric HisF-based 
TIM barrel!?°l371 have also been created. 


The existence of 4/8-fold internal symmetry was sug- 
gested by Soding et al."3! based on a computational 
analysis of TIM barrel sequences. For example, Esche- 
richia coli KDPG aldolase!*! (PDB: 1FQ0) was suggested 
to possess a distinct 4-fold symmetry, with discernible 
8-fold symmetry. The design of a 4-fold symmetric TIM 
barrel"! confirmed the possibility of higher orders of in- 
ternal symmetry in natural TIM barrels, and will be dis- 
cussed in detail in the next section. It should be noted 
that no experimental evidence for the existence of 8- 
fold symmetric TIM barrels has been reported to date. 


De novo TIM barrel design 


The TIM barrel fold has been a long-standing target for 
de novo protein designers. As previously described, nu- 
merous TIM barrels have been successfully designed 
based on preexisting natural half-barrels. In contrast, 
the de novo design of TIM barrels occurred in incremen- 
tal steps over a period of 28 years.!**! 


The Octarellin series!o1471481149150] of proteins (Octarel- 
lin |> VI) were the first attempts to create a de novo TIM 
barrel. As the field of protein design was still in its in- 
fancy, these design attempts were only met with lim- 
ited success. Although they displayed circular dichro- 
ism spectra consistent with aB proteins and some coop- 
erative folding characteristics, all Octarellin series pep- 
tides were insoluble, and had to be resolubilized from 
inclusion bodies for further characterization. Interest- 
ingly, Octarellin V.1°! displayed a Rossmann-like fold 
under co-crystal conditions. 


The Symmetrin series!” of proteins (Symmetrin-1—4) 
displayed more favorable biophysical characteristics. 
Symmetrin-1 was readily soluble, displayed circular di- 
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Figure 8 | sTIM-11, the first successful de novo TIM barrel 
design. The asymmetric (aB)2 units are colored distinctly, 
highlighting the internal 4-fold symmetry. 


chroism spectra consistent with aB proteins, and dis- 
played excellent cooperative unfolding and refolding 
characteristics. Despite these advances, all proteins in 
this family displayed molten characteristics when ana- 
lyzed using NMR (nuclear magnetic resonance), and 
further work to solve their structures could not be pur- 
sued. 


Proteins of the sTIM series’! represented the first suc- 
cessful de novo TIM barrel design. sTIM-11 (PDB: 
SBVL) was designed with an internal 4-fold symmetry, 
to reduce the complexity of computational design using 
the Rosetta software suite.?! Previously-derived first 
principles“! were used to delineate secondary struc- 
ture topologies and lengths. sTIM-11 proved to be a 
highly thermostable, cooperatively folding design that 
adopted its intended structure (Figure 6). 
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